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We microscopicaUy derive transport equations for the conduction electrons in ferromagnetic ma- 
terials with an inhomogeneous magnetization profile. Our quantum kinetic approach includes elastic 
scattering and anisotropic spin-flip scattering at magnetic impurities. In the diffusive limit, we cal- 
culate the resistance through a domain wall and find that the domain- wall resistance can be positive 
or negative. In the limit of long domain walls we derive analytical expressions and compare them 
with existing works, which used less general models or different theoretical frameworks. 
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I. INTRODUCTION 

Conducting magnetic materials are an active research 
topic at present due to promising applications like mag- 
netic memory storage devices which make use of magneti- 
zation reversal in pillar multilayer nanostructures^i^i^i^ 
or domain wall motion^'^'^'^ i^^'^"'^ as proposed for the race- 
track memory-'^^. On one hand, domain wall motion is 
realized by sending spin-polarized current through the 
domain wall, so that the mutual interaction of the elec- 
tron spin with the ferromagnetic order parameter leads 
to a motion of the wall. This is due to the so called spin- 
torquei^ii^ii^ii^iii, the transfer of spin-angular momen- 
tum. On the other hand, the electronic current flow is 
also affected by the presence of an inhomogeneous mag- 
netization. Most prominently, there is a change in the 
resistance when the current runs through a domain wall 
in comparison to the resistance in the absence of the do- 
main wall. The resistance change can have different ori- 
gins that can be seperated into the extrinsic and intrin- 
sic domain- wall resistance (DWR). The former includes 
orbital and anisotropic magneto-resistance. The latter 
contains the direct influence the domain wall has on the 
electronic conduction channels: if the magnetization di- 
rection is not homogeneous in space, the spin majority 
and minority channels are no longer eigenstates, which 
in turn changes the conduction properties and also can 
have influence on the impurity scattering rates. There is 
also spin accumulation in the vicinity of the domain wall 
which leads to an additional potential drop. In any case, 
the extrinsic mechanisms have to be carefully identified 
in order to obtain the intrisic domain- wall resistance from 
experiment. The DWR has been studied in a number of 
works in the past, both theoreticallyi2iiaiSai2ii^i2^iSii2^i2£ 
and experimentally^li^i^Sii^S. Reviews about DWR in 
nanowires made from ferromagnetic transition metals, 
experimental measurements and details on the treatment 
of extrinsic magneto resistance can be found in^LiH. 

On the theoretical side, several limiting case have been 
investigated using a variety of theoretical methods. The 
worksi2i2£iiSi perform a diagrammatic evaluation of the 
Kubo-formula introducing scattering in the unperturbed 



Greens functions by two phenomenological parameters 
T^l , the momentum scattering times for spin up and down 
channels. In this calculation, spin-flip processes are not 
included. As we will discuss later, this leads to a spin 
accumulation that does not decay even arbitrarily far 
from the domain wall. Hence, this neglect of spin-flip 
is only possible, if the distance between the domain wall 
and leads is much smaller than the spin-diffusion length. 
In a complementary approach. Levy and Zhangi^ use a 
linearized Boltzmann equation. They do not consider 
changes in the electronic spectrum, i.e. they assume spin- 
independence of the wave vector = fcj , restricting the 
validity to the regime of small exchange splitting. Their 
analytical calculation is done in a basis that diagonalizes 
the Hamiltonian, which is possible in case of a constant 
magnetization gradient, known as spin-spiral. Thus, they 
cannot take into account a finite contact geometry and 
finite domain wall length, but have to consider an in- 
finitely extended spin-spiral for which they calculate the 
conductivity. Spin-flip processes are absent, so, again, 
the above statement concerning spin accumulation ap- 
plies. Furthermore, they perform a multi-pole expan- 
sion of the distribution function but only include terms 
up to the p-wave component. However, as we will see 
during our calculation, this is not sufficient in general. 
Lastly, we believe the Boltzmann equation, they use, 
lacks terms that should appear as a result of the gauge 
transformation. Bergeret et al use the Keldysh technique 
to derive a quasiclassical equation valid in the diffusive 
limit^. However, they consider a different regime of va- 
lidity, in which the scattering mean free path 1^ is the 
smallest length scale in the system (besides the Fermi- 
length) , and not the precession length as will be the case 
in our treatment. Likewise, they do not consider spin- 
flip processes, even though during their calculation, they 
perform steps which implicitly require longitudinal spin 
excitations to relax. Finally, Simanek et al?^'"^^ used 
equations of motion for the quantum distribution func- 
tion in Wigner space, which however contain a term that 
we cannot reproduce. Before, Bergeret et al. noted that 
this term violates particle conservations!. Nevertheless, 
this term does not affect the statement of Simanek et al. 
that there is quenching of the spin-accumulation due to 
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rapid transverse precession. This also emerges from our 
theory and we will make use of it later (see the discussion 
around Eq. ([70| ). 

In this article, we pursue a fully microscopical the- 
oretical approach to the DWR in the limit of wide 
walls, so that quantum mechanical electron reflection 
at the domain wall can be neglected. This allows us 
to use a standard quasiclassical approximation and ne- 
glect spin-dependent scattering due to abrupt potential 
changes^^^i^. In section II, we begin by introducing our 
model and deriving a quantum transport equation using 
the Keldysh kinetic equation approach. These provide a 
rather general framework to investigate a large variety of 
static transport problems. In Section III we solve these 
resulting equations analytically in certain limiting cases 
for model domain walls and discuss our results and re- 
late them to various existing theoretical works dealing 
with the issue of DWm^^^^^^^^^^^^^^. Finally we 
conclude with an outlook on open problems. 



II. QUANTUM TRANSPORT EQUATION FOR 
FERROMAGNETIC CONDUCTORS 

In this section we derive a quantum transport equa- 
tion from a model Hamiltonian that describes the kinetics 
of conduction electrons in materials with inhomogcncous 
magnetization profile. 



A. Model and Hamiltonian 

We consider a system of effectively non-interacting 
electrons whose spin degrees of freedom are coupled to 
the ferromagnetic order parameter in the mean-field ap- 
proximation via the spatially dependent exchange field. 
The single particle Hamiltonian has three contributions, 

n = Ho + ns + v (1) 

= {^j^-SfJ-ir)~etp{r)ji-—m{r)& + Vir) 

Here, Tip is the usual free quasi-particle energy contribu- 
tion with the dispersion relation Ck = 3^ and effective 
mass M and in spatial representation, k = — zV^- ^{r) is 
the external electric potential felt by the quasi-particles 
of charge — e. 6^ is a. shift in the chemical potential due 
to the magnetization gradient drtn which later turns out 
to be of order (5/i — 0{drm)^. TCs describes the coupling 
of the electron spin to the exchange field with a constant 
magnitude A and the local magnetization direction de- 
noted by the unit vector m{r). The 2x2 matrix spin 
structure is denoted by a- = (ai, a'2, 173) is the vector 
of Pauli matrices, such that the electron spin operator is 
given hy S = ^Sr. 

In accordance with the mean-field approach, the length 
scale of the spatial variations is much slower than the 



relevant atomic scales. More specifically, this condition 
reads 

l/lm = \drTn\ ^Pf/H . (2) 

The exchange field is created by electrons that align their 
spin preferably in the same direction due to the (here 
ferromagnetic) exchange interaction. In conducting fer- 
romagnets, the electrons contributing to the local mag- 
netization can either be localized and, thus, do not par- 
ticipate in transport (d-electron character) or be delocal- 
ized and, hence, are subject to electronic transport phe- 
nomena (dominant s-electron character). These extreme 
cases constitute two distinct models with the major dif- 
ference being the way in which the self-consistency condi- 
tion for the exchange field is employed. These are known 
as s-d model and itinerant Stoner model, the latter one 
describing a system in which transport and magnetism 
arise in fact both from the same delocalized electrons. 
However, real physical systems are usually between these 
two cases. Below, we will restrict ourselves to the s-d 
model in which the magnetization profile remains static 
even if the conduction electrons are in a non-equilibrium 
configuration. Note, that fluctuations of the order pa- 
rameter are neglected. 

We disregard the influence of the effective magnetic ex- 
change field on the electronic orbits, which represents the 
Lorentz force and leads to the orbital magneto-resistance 
(OMR). Theoretically, as well as experimentally, the 
OMR and other extrinsic contributions such as the AMR 
(anisotropic magneto- resistance), which stems from spin- 
orbit coupling and leads to a resistance depending on the 
angle between magnetization and current directions, can 
be separated from the true DWR and therefore are not 
included in this article. For Bloch walls and the CPW- 
geometry (current perpendicular to wall), which we em- 
ploy in our model later on (see FIG. [T|), the AMR plays 
no role, since the magnetization direction is always per- 
pendicular to the current direction. For other walls like 
Neel walls, the spin direction within the wall attains a 
component parallel to the current and the AMR has to 
be carefully distinguished from the DWR. 

The impurity scattering potential can be divided into 
two contributions, 

V(r) = l/i(r)-HKnag(r) . (3) 

Vi describes the scattering from randomly distributed 
static impurities and for point-like scatterers has the 
property 

(M(r)M(r'))imp = X^'^'5(r-r') , (4) 

equivalent to the treatment of 14 as a delta-correlated 
fluctuating Gaussian field, x'"^^ measures the strength of 
the impurity scattering potential and ()imp denotes the 
averaging over all impurity configurations. 

In a similar way, scattering at impurities that have in- 
ternal spin-degrees of freedom makes the scattering ver- 
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tex spin-dependent, such that 
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(V^mag(r-)o\/„iag(r-'))imp = Yl xi]''^ ° ^j^ir - v') . (5) 

Therefore, magnetic impurity scattering in the point-Uke 
hmit can be treated as delta-correlated fluctuating Gaus- 
sian magnetic field which couples to the electron spin by 
the usual Zeeman term. The size of the fluctuations can 
be generally spin-anisotropic which manifests itself in the 



tensor structure of x. 



(m) 



B. Keldysh Technique 



Greens function. For the remaining part of this article, 
we set h — 1. 

The further treatment is done in the Wigner- 
representation which is obtained from usual spatio- 
temporal representation via the transformation 



G{k,e,r,t) = / d'^zdr e 



-ikz+ieT 



xG(r + -,^+-;r--,t-- 



(6) 



The standard way to proceed in non-equilibrium 
physics is to set up the kinetic equation for the Keldysh 



The product of operators has to be carried out using the 
rule 
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C{k,e,r,t) = A{k,e,r,t) ® B{k,e,r,t) = e, r, t)e^(^'- ^'""^'^^ 

I 



? k 3t-+ d edt) 



B{k,e,r,t) 



(7) 



where d and d denote derivatives acting only to the left 
and right, respectively. The transformation ([6]) intro- 
duces center of mass coordinates r, t and Fourier trans- 
formed relative coordinates k and e, respectively. Note, 
that the product (g) is associative. 

The Greens function G is defined as expectation value 
of the electron field operators, time-ordered along the 
Keldysh contour'^^. The ordering along backward and 
forward time Keldysh contour gives rise to an additional 
2x2 matrix structure (denoted by " ) , which in an ap- 
propriate basis takes the convenient form 



G^ G< 
G^ 



(8) 



and G^ are the retarded and advanced Greens func- 
tions, well known from equilibrium theory and carry in- 
formation about the spectrum of the system. In particu- 
lar, one obtains the spectral density simply from 



R 



G^ 



(9) 



The spectral function generally obeys the normalization 
condition 



— A{k,e,r)^l. 
Zn 



(10) 



The fc-integration of the spectral function yields the den- 
sity of states, here defined as number of states per unit 
energy and volume: 



(27r)3 



A{k, e, r) — 2ni'{e, r) 



(11) 



The lesser component G^{k,e,r,t) describes the occu- 
pation of states of the many particle system and is given 
in terms of electron field operators by 

G<^(r,r',t,t') = *(v4(r',t')V'a(r,t)) , (12) 

where the grandcanonical average is taken and the indices 
a, (3 are one of t, i- In equilibrium, it takes the form 

G<(fc,e,r) = zi(fc,e,r)/D(e) (13) 

with the Fermi-Dirac distribution function 

1 



./d(6) 



(14) 



where f3 — denotes the inverse temperature. 

From the lesser Greens function G^ we can easily ob- 
tain various physical quantities of interest such as the 
quasi-particle spin-charge density 



N{v,t) 



(27r)3 2^1 
and spin-charge current density 



— G<(fc,e,r,i) , (15) 



J(r,i) = -e 



(2^y_2^'^^^'(^'^'^'^)' (1^) 



where the quasi-particle velocity is Vk ~ dk^k = k/m. 
The spin-charge density iV is a 2 x 2 matrix in spin space 
and its trace yields the charge-density, ric = TriV while 
the spin-density is given by s = — ^Tr(<TiV). In a similar 
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manner, we obtain the charge current = Tr j and the 
spin-k current j^. — — :|jTr((Tfc J). 

In order to find G for a specific physical system, we 
need an equation of motion, called the Dyson equation, 
which can be written in the two forms 



(ei -Ho-ns)i-t 



0(5 = i , 



Gi 



(ei - Tio - Hs) i - S = i 



(17) 



The self-energy S incorporates scattering by magnetic 
and non-magnetic impurities. The leading contribution 
is calculated in the self-consistent Born approximation, 
which truncates the series of irreducible diagrams due to 
multiple impurity scattering after the first one: 



(VGV) 



imp 



(18) 



In physical terms, this means that all kinds of interfer- 
ence effects like weak localization are dropped from the 
theory. In this approximation, the self-energy for spin- 
independent impurity scattering takes the form 



Si 



(27r)3 



G. 



(19) 



Similarly, for magnetic impurity scattering we obtain 



^mag 



(2^ 



(T i G(7 j 



(20) 



The tensorial structure of ^(™) accounts for situations in 
which scattering is anisotropic in spin space. For exam- 
ple, if an exchange coupling between the internal impu- 
rity spin and the ferromagnetic order parameter exists, 
the spin is preferably aligned along this direction. Con- 
sequently, the impurity will scatter the electron with a 
magnitude depending on its spin. In the case of uniaxial 
symmetry (the symmetry axis is denoted by n). 



(m) 



xi"^(4 



^ I ("') 



(21) 



For the above example, the unit vector n actually cor- 
responds to the local magnetization direction m. We 
assume throughout the rest of this article that scattering 
is weak such that 



(22) 



where is the Fermi energy of the conduction electrons. 
Put in other words, this means that transport quantities 
like density of states do vary very slowly on the scale of 
the relaxation rates given by S. This assumption will 
allow us later to make use of the quasi-particle approxi- 
mation. 

C. Spectral Density 

The first step in solving the Dyson equation (fT7|) is to 
find a solution for the spectral function. To do so, we 
write equations for the retarded/advanced components 
of G, 

[e-€k + 5^l{r)]i + ^m{r)&^ ® G^'^ = I , 
G^'^ ® (^e -€k + 5^l{r)]i + ^m{r)a^ = 1 . (23) 

(^R/A only differs by the boundary condition which can 
be simply incorporated by substituting e e ± 177 in 
Eqn. (|23|) . Furthermore, since S <Si E-p, we will not in- 
clude impurity self energy corrections to the spectrum, 
i.e. we neglect the line broadening and assume a delta- 
peaked spectrum, commonly referred to as the quasi- 
particle approximation. Also, the external electric field 
affects the spectrum only in second order of the field. 
To proceed, we perform a gradient expansion 



l + -idrdk - dtde 



dkdr + dedt 



(24) 



and determine G^/^ iteratively by the order of the gra- 
dient. It turns out, that we need only up to order d^. For 
example, one finds for the first two terms of the spectral 
density A = Aq + Ai + 



Ai 



= 2tt 



& {mx{vkdr)m) [de S{e - el) + S{e - ej,) 



2 
A 



5{e 



(25) 
(26) 



Here, we defined the projectors in spin-space P|^| = ^(i± 
ma) that project on the spin up/down direction and 
e'y} = efc =F 3- is the dispersion relation for majority and 
minority spin bands which are exchange split due to the 
s-d coupling. 



The density of states for majority and minority spin 
bands are 



ie) 



d^k 
(27r)3 



6{e 



(27) 
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With i^i i without energy argument we denote the density 
of states at the Fermi-level. The matrix density of states 
without magnetization gradient present is thus 



V = \{l + ma)v^ + \{i-m&)vi 



va{i+Pm&) , (28) 



which leads to 



(30) 



where we have introduced the polarization of the Fermi 
surface 



P 



(29) 



and vq = \{v] denotes the average density of states. 

The term in Eqs. (j23p constitutes an additional 
chemical potential shift due to the magnetization gra- 
dient and is determined by enforcing local charge neu- 
trality. The total quasi-particle density is constant 
throughout the ferromagnet, i.e. Tr7V(r) — eiVion, 



J 



D. Kinetic Equation for the Non-Equilibrium 
Distribution 



To obtain a kinetic equation for the non-equilibrium 
distribution function, we subtract the left- and right- 
conjugated Dyson equations (fT7|l . 



The equation for the lesser component is 



(31) 



(e - Cfc + ^A* 



\{A-t<]-\{v-G<] 



r 



+ i 



KG? 



(32) 



where the (anti-)commutators are defined by {A ® _B} — 
A® B + 13 ® A and [A® 13] = i B - B i. The 
imaginary part of the self-energy 



(33) 



describes relaxation due to impurity scattering, while the 
real part is responsible for changes in the energy disper- 
sion relation. There exists a Kramers-Kronig relation 
between real and imaginary part. 



1 



3?S(e) = -(E 



ZTT 



de' 



, r(e') 



(34) 



However, as can be seen from the integral representa- 
tion in this formula, the real part depends on the com- 
plete electronic spectrum of the system, since the scat- 
tering rate is directly related to the density of states, 
r(e) oc v{e). Considering that the dynamics accompa- 
nied by a rotation of the magnetization direction, be it in 
time or space, is affecting only an energy region of the or- 
der of A, these changes constitute only a tiny fraction of 
the whole energy range. Thus, corrections due to magne- 
tization dynamics to the real part can be neglected when 
compared to the whole background contribution, which 
then is just a constant (however formally diverging due 
to the assumption of fe-independent impurity scattering) 
and merely renormalizes the electronic spectrum. In fact, 
the only reason to include impurity scattering is to add 
momentum and spin relaxation to the conduction elec- 
tron system, as contained in the imaginary part of the 
self-energy, F. Therefore, the two commutators involv- 
ing real parts can be dropped from equation 



Generally, we can distinguish the contributions to the 
kinetic equations emerging from two regions in energy. 
The electronic states deep inside the Fermi sea are af- 
fected by an inhomogeneous exchange splitting Am(r). 
The fraction of the Fermi sea that contributes is given 
by A/i?F and not necessarily small. However, due to 
the Pauli principle, these states are fully occupied for 
all reasonable temperatures and the change in the spec- 
trum docs not affect the dynamics of the mobile elec- 
trons close to the Fermi surface. These are responsible 
for the dynamics in the quantum kinetic equation, since 
the electrons in an energy window given by the temper- 
ature, voltage or other low-energy scales have the free- 
dom to move. These differences can be used to eliminate 
the high-energy contribution from our quantum kinetic 
equation. We perform a series of three steps to obtain an 
equation describing the low-energy dynamics alone. 

First, we eliminate the electric potential by the sub- 
stitution u> = € -\- eLp, which transforms the derivatives 
according to 



, dr ^ dr ~ eEda. 

with the electric field E — —drf- 
Secondly, we make the ansatz 



(35) 



G<(fc,w,r) = iA{k,uj,r) fo{uj-eip) + SG{k,uj,r). (36) 



The first term drops out from the kinetic equation 
leaving an equation for 6G alone. SG has two very prac- 
tical properties. It is proportional to the electric field, 
which allows us to drop the term eip{r)dG, since we are 
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interested only in linear response to an external field. 
The external potential is incorporated via appropriate 
boundary conditions that are concretized below. Fur- 
thermore, 5G is peaked around the Fermi-level, reflecting 
the fact that non-equilibrium processes take place only 
in the vicinity of E-p, provided the temperature is low 
enough. In fact, we use the zero-temperature approxi- 



3/d 



-5{e — E-p) throughout this work. 



mation 

This brings us directly to the third step which consists 
of integrating the whole equation over energy after set- 
ting Lo = Ep in all prefactors to 5G on the right-hand-side 
of equation ([5^ . The spectral densities in the collision 
integral become A(fc, Ep, r), which means we neglect the 
energy-dependence of the scattering rates. This is per- 
fectly compatible with the linear response regime, since 
a more in-depth investigation shows that corrections due 
to the energy dependence of A are of quadratic order in 
the electric field. 

For the stationary situation {pt = 0), the resulting 
kinetic equation for 



g{k,r) 



(37) 



finally reads 
A 



Vkdrg - 



I— [ma- ° g] + {drSn)dkg = Xi[g] +T^[g] ■ (38) 



The collision integral takes the form 



^[A{k,Ep,r)° g{k',r)} 



(27r)3 

i{i(fe',i?F,r)° g{k,r)} 



(39) 



and similar for the magnetic scattering contribution 
2rn [.9] • ° denotes the Wigner product ([7]) with derivatives 
dr and dk only, and does not act on k' . Again, we will 
need to include only terms up to order . Equation ([55]) 
along with (|39p constitutes a linear integro-differential 
equation for g and serves as the basis for our further an- 
alytical treatment. 

Provided we have found a solution for g with appro- 
priate boundary conditions, we can finally determine the 
physical observables of interest by simply substituting 
the ansatz into In this way, we obtain the 

quasi-particle spin-charge density 



/do; f d^k 1 f d'^k 

2^ J (2^^^^"' Mi^)-eMr)^ J j:^ A{k , Ep , r) + n{r) (40) 

I 



with the low-energy spin-charge density excitations n = 
— e / (27r)^a ff- We made use of the zero-temperature ap- 
proximation and smallness of the external potential in the 
linear regime, viz fpi{uj — eip) = fpi{uj) + eip 5{ijJ — Ep). 
The (spin-)current is obtained in a similar straightfor- 
ward manner. 

On length scales much larger than the Thomas-Fermi 
screening length, which is of the order of atomic distances 
in metallic materials, local charge neutrality is fulfilled. 
Since the first term on the right-hand side of ([^0]) corre- 
sponds to the equilibrium density which exactly neutral- 
izes the material, the trace of the last two terms has to 
vanish, 

Trfi(r) = e\{r)^ f ^Tri(fc, Sp, r) . (41) 

Therefore, Trn(r) is directly related to the local electric 
potential and thus the local electric field E{r) — ~dr^{r) 
in the system. 

Now we consider a region far from any inhomogenities 
in the magnetization for which we want to specify an ap- 
propriate boundary condition. In the absence of a mag- 



netization gradient and non-equilibrium spin-excitations 
in the conduction electron system, a less strict equality 
between the last two terms of (|40p holds, 



n(r) = eV(r)^ J -^^Ml^, Ep, r) . (42) 



In a simple onc-dimensional geometry and at the left and 
right boundaries xl.r, the boundary condition simply be- 
comes 



ni^L,R) = eV(a^L3) (i'tA (43) 



Before continuing to solve the kinetic equation for 
g, let us specify a convenient geometry for the contact. 
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FIG. 1: Contact of length L with a Bloch domain wall of 
characteristic size w situated in between. The domain wall is 
described by the magnetization gradient qix) = dQix). The 
lower plot shows my{x) = ~q{x) and mz{x) according to 
equation (|47p with unity wall width w = 1. 



III. DOMAIN- WALL RESISTANCE IN 
DIFFUSIVE WIRES WITH CURRENT 
PERPENDICULAR TO WALL (CPW) 
GEOMETRY 

A. Model of the Contact 

We treat a ferromagnet in a quasi one-dimensional ge- 
ometry such that there are only gradients of the magneti- 
zation in the x direction, while in the y- and z-directions 
the system is homogeneous. Furthermore, we assume a 
coplanar magnetization which allows to parameterize the 
magnetization direction by a single angle 9{x), defined by 

( ' \ 

m{x) = sin6'(a;) . (44) 

\ cos0{x) J 

Current flow perpendicular to the domain wall implies 
that the direction of j is along the z— axis, so that the 
problem becomes effectively one-dimensional. Therefore, 
in this case, every quantity of interest only depends on 
X, we can set dy = = and abbreviate d = d^- 

The contact of length L contains a domain wall of char- 
acteristic width w. The left- and rightmost parts have 
opposite magnetization directions and a so far arbitrary 
magnetization profile 9{x) connects the two magnetic do- 
mains. The situation is sketched in Fig. [TJ 

It is then convenient to perform a local SU(2) gauge 
transformation in order to rotate the coordinate system 
used to represent the spin-orientation such that its z-axis 
is always aligned to the local magnetization direction m. 
The unitary operator [W = IA~^) 

U = e-^e*- (45) 

exactly performs this rotation in spin-space using the pre- 
viously introduced angle 0, so that U mcr = ^z- As 



9 is spatially dependent, the derivatives obtain an addi- 
tional term after the transformation, 

U d U-^ ^ Id + ^q{x)cr^ = d . (46) 

Henceforth, we have to deal with the gradient q{x) = 89 
in our equations and which fully parameterizes the do- 
main wall. The gauge transformation (|45p rotates the 
orthonormal basis system {v xm^v ,m) such that it be- 
comes the canonical set basis vectors, {ex,ey,ez)- The 
unit vector v points in the direction of change of the 
magnetization m, i.e. dm — qv. Note that v is perpen- 
dicular to m, since m is a unit vector. 

An analytical form of the domain wall profile can be 
derived in the context of minimising the free energy of 
the ferromagnetic material, commonly a compromise be- 
tween exchange and anisotropy energy. Typically, one 
finds 

cos6l(x) = tanh(^) , (47) 

so that 

q^x) = ^ = tanh(i) , (48) 

where w constitutes a typical length scale over which the 
domain wall extends. 



B. A Hierarchy of Equations 

In light of the solution procedure that is to come, it 
is convenient to use a 4-component vector representation 
(which is denoted by an arrow to distringuish it from the 
3-component real-space vectors set in boldface) , in which 
the spin-charge density excitations and current density 
take the form, respectively, 

n = (n+,n_,n|,n|) , 
J = U+J-Jhh) ■ 

Starting from the 2x2-matrix reprentation after the gauge 
transformation, h = ^l + n^ax -\-nyay + nzd'z, we define 
the spin-up/down densities — ^ iriz- Additionally, 
the transverse spin-degrees of freedom are transformed 
according to n± = ± iuy, which corresponds to circu- 
larly polarized transverse spin excitations. This basis is 
convenient because it diagonalizes the equations of mo- 
tion for the transverse dynamics in absence of magneti- 
zation gradients, see Eqn. ((56|) and following. 

The full kinetic equation for g is still very involved, so 
we need to perform additional simplifying steps. There- 
fore, we multiply Eq. ([55)) with {n > 0) and afterwards 
integrate over the whole fc-space. This involves evaluat- 
ing terms of the form 
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and defining moments of the Greens function g 



(49) 



The first two momenta are 5*^"^ — n and g^^-* = j, ob- 
viously. The kinetic equation psp turns into an infinite 
hierarchy of equations relating these moments. It is nev- 
ertheless possible to find an analytical solution to these 
equations in form of a systematic series expansion in the 
function q{x). As we will show below, it is then possible 
to truncate the hierarchy by setting g^^^ = in order 
to calculate the domain-wall resistance. The first two of 
these equations are reminiscent of spin-charge continuity 
and diffusion equations and take the explicit form 



(50) 
(51) 

where we defined the derivative d, modified by the SU(2) 



Tn + d^i = 0{q) 
n[j + D^dri\ = 0{q) 



gauge transformation (see equation 



d=ld + = ld+-q 



/ 1-1 

0-11 
1-10 

V-1 1 



(52) 



The right-hand-sides of these equations include impor- 
tant contributions that are due to modification of various 
transport properties in presence of a magnetization gra- 
dient like the change in the density of states which enters 
the collision integrals. Thus, the right-hand-side vanishes 
as g — > 0. Their explicit form along with the equations 
for g^^\ etc. are rather lengthy and not needed for our 
following discussion. Further details can be found in the 
appendix in equations (|A17[) - (|A25[) . 

The various scattering rates associated with momen- 
tum relaxation H and spin-flip/-dephasing F are most 
conveniently expressed, introducing the parameters 7, k 
and ^, in the following way: 



n = diag (7,/- + iA, 7^ - zA,7 



7. 


V 















= Tj^diag ( 1 + -, 1 - -, 1 + 7, 1 - 7 






In In 



I 

X 

i 



X 






1(1 -P)^^ 



\ 



-^{1 + P)k 
i(l + P)« / 



Spin precession is incorporated as well and manifests it- 
self as an imaginary part in the entries of the transverse 
subspace of F and H. 

In terms of the parameters of our specific model, the 
various scattering parameters take the form 



^] X^' +X\ + 2x1 



1 = P 



X«+x[r^+2xf^ 
2xi"' 



(m) (rn) 

X I +X1 



(53) 



X measures the ratio of momentum relaxation rate to 
exchange splitting A, ^ denotes the transverse spin- 
dephasing strength and 7 can be considered as scattering 
asymmetry parameter while k, is the strength of spin-flip 
scattering. 

Of course, the parameter 7 carries an implicit depen- 
dency on the exchange splitting A entering through the 
polarization parameter 



P(A) 



1 + Vl - A2 ' 



(54) 



expressed here by the dimensionless exchange splitting 
A = 2^ ■ Within our specific model of impurity scatter- 
ing, 7 — P(l — ^k) is not independent of the other param- 
eters. It is worth mentioning that this specific form of 
7 has to be employed in order to guarantee for a consis- 
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tent treatment of the gradient corrections in the coUision 
integral. 

Furthermore, the matrix of diffusion constants reads 



R = diag(— p-,;;^ :,D^,Di 



-D_Ldiag 




l + A 1-A 



X + iX — i 1+7 I-7 



(55) 



where we introduced Dj_ — fff-^ — ^''^f^- Transverse 

spin excitations are also subject to precession as they 
diffuse which results in the complex values of the effective 
transverse diffusion constants. 

Finally, let us stress that the reason for the inclusion of 
the moments up to g^'^^ lies in the fact that in our regime 
of investigation the precession length Zp^c = "^^^/sa 
besides the Fermi-length, the smallest length scale in the 
system. As can be seen from the equation (l60l) . Zproc 
describes the period of oscillation of transverse spin ex- 
citations. In the diffusive approximation, as for example 
used by Bergeret et aP4, the scattering mean free path 
Is is the smallest length scale in the system (besides the 
Fermi wave length) and not Ip^ec- Hence, it is possible to 
truncate the hierarchy of equations already by 5^^-' = 0, 
which results in only two equations that are just the spin- 
chage continuity and diffusion equation. The latter is ob- 
tained by plugging the equation for g^^"^ into the equation 
for g^^' (see also appendix, for more details). 



straightforward manner from equations (|A17P - (jA2ip by 
collecting all terms associated with the corresponding 
factor d^qd^qd''. Restricted to terms that contribute to 
DWR, its explicit form (see eqn. (|A23[) ) is given in the 
appendix. 

Since we have a perturbative treatment in q in mind, 
we determine the Greens function of Eq. (I56p 



(r - Dd^) G{x) = lS{x) 



(58) 



Separated into longitudinal (1) and transversal (t) sub- 
space, the Greens function is 



G(x) 

Gii(x) 
Gtt(x) 



Gtt(x) 
Gu{x) 



He 











X 



(59) 
(60) 



X 



e = 



' 2n 



1 J_p^k\^\ 



2k 

2 



prcc 



(1 - ix) (1 - txO 



Matrices in the 2x2 subspaces are denoted by a sin- 
gle underbar as compared to the double underbar which 
indicates a 4 x 4 matrix. The index 1 and t refers to lon- 
gitudinal and transverse componenents, respectively, so 
that for example 



C. Solving the Hierarchy of Equations 



We eliminate the higher order moments by iter- 

atively substituting the equations into each other, care- 
fully keeping terms that contribute up to order q^. We 
find a differential equation for the vector of quasi-particle 
excitations n of the form 



(r-Dd^)n 



Wiq)n 



(56) 



where the differential operator W_{q) vanishes for g — > 
and contains all possible terms up to order q^. In the 
homogeneous {q = 0) and coUinear case, Eq. ([55)1 cor- 
responds to the transport equation used by Valet and 
Fert'^'^-. Let us stress that n represents the unscreened 
spin-charge density excitations, while the true, screened 
quantity n'-^^ is, according to Eq. (|^ . 



(57) 



Here, spin-excitations are not screened since our model 
does not include a spin-dependent interaction. 

Generally, WJ^q) contains terms of the form 
Yjjkd^qd-' qd'^ , Y^ijd^qd^ and Y^id^, where the order 
of q and d is crucial, since d acts on everything to 
its right. The constant matrices Y_ that depend on 
our set of parameters, A,7,^, k can be obtained in a 



The longitudinal component G\\ consists of two contri- 
butions. The first term of Gji decribes spatial damping 
of spin-up/down non-equilibrium excitations which man- 
ifests itself in the characteristic exponential decay on the 
spin-diffusion length. 



1 

A2 



2A. 



7i 



(61) 



The second term of G\\ describes the linear behaviour of 
the chemical potential in a homogenous system in the 
absence of a magnetization gradient. The two tensors 



H 



1 



K = - 



2 {DtD 
1 



Dili 



obey the useful identities 



In in 
In In 



TuK 



1 



Fn - 



H = 



A 



DniR^K) = 1 



(62) 
(63) 

(64) 
(65) 

(66) 
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which can be invoked to easily verify that the longitudi- 
nal Greens function G\\{x) in fact fulfills equation ([55)1 . 
The real part of the complex wave- vector k describes the 
precession of transverse non-equilibrium spin excitations 
while its imaginary part is the damping due to dephasing 
mechanisms. Therefore, the root of k has to be chosen 
such that is has a positive imaginary part. 

In the following, we restrict ourselves to the regime 
X <C 1, viz. a momentum relaxation rate much smaller 
than the exchange splitting. Again, since leading order 
correction term turns out to be of order , we drop any 
terms of higher order than that. Later, we will see that 
this restricts the validity of the result for the DWR to do- 
main wall lengths larger than the spin-precession length 

^prcc ■ 

In this regime, the transverse oscillations are very rapid 
on the scale of the magnetization gradient. Therefore, it 
is suitable to eliminate the transverse degrees of freedom 
by first splitting the equation of motion (|56p for n into 
transverse and longitudinal parts, 

{r-Dd^)^^nt = Wtm + WtSt (67) 
(E-^c*')!!"! = Wnni + Wurit (68) 

and writing down the formal solution for the transverse 
component 

/oo 
da;'Gtt(a; - x') [Wtini{x') + WttUtix')] ■ 
-OO 

(69) 

In the limit x ^ li Gtt{x') varies on a scale determined 
by ^ = Zprcc which is much smaller than other length 
scales of interest, which are variation of m and the ex- 
ternal electric field. In particular, for the length of the 
domain wall, w 3> Zpioc- Thus, to leading order in x> 
we can consider Gtt as a representation of the Dirac S- 
function and perform the integration. We obtain 

nt{x) = Z[I£ti"i(a;) +l£tt^t(a;)] . (70) 

Here we introduced the spatially integrated transverse 
Greens function, 

dxGu{x)=(^^ ) , (71) 

where 

F - « -^^X^ = -4 ■ (72) 

Note that £ simply corresponds to the inverse of the 
transverse part of F, a fact which becomes clear by noting 
that our approximation corresponds to the neglect of the 
transverse diffusion term. 

Additionally, we only need to keep the first term of 
(I70p since the backaction on the transverse dynamics, 
represented by the second term, appears only in higher 



orders in q and x- Explicitly, this is expressed by the fact 
that to leading order in q, W_tt vanishes, so that 

nt{x) = EWt\n\{x)+EWttEWt\n\{x) + ... 

= E WtA{x) + O(x') (73) 

Putting this result back into the equation for the longi- 
tudinal dynamics yields the formal solution, 

/C30 
dx' Gn{x' - x) (74) 
-OO 

X [mini{x') +mtZWtMx')] . 

The boundary condition (|43| for the left and right side 
of the contact reads 

n(±f) = eV(±f) (0,0,i^T>^i) ■ (75) 

The zeroth order solution with externally applied bias 
voltage V , that satisfies this boundary condition, is sim- 
ply found to be 

n'^^){x)^^^Ex (0,0,1^1,^.4) , (76) 

where E — \s the constant external electric field in 
absence of the wall. 

Substituting 'n!^'^\x) into the right-hand-side of the 
solution ([74|) yields the second order correction in g, 
n^'^\x). Due to charge-current conservation, the cur- 
rent flowing through the contact is still determined by 
because 9n'^'(x) is taken to vanish at the bound- 
aries by assuming that the contact is long enough for 
any finite size effects to become negligible, i.e. L ;2> A. 
In this regime, the exponential term in the longitudinal 
Greens function ([59]) can be neglected which means that 
spin-accumulation has faded near the reservoirs. Then, 
the current can be deduced directly from Eqs. (|76p . un- 
affected by the correction n'^^ Explicitly, this current 
reads 

J(o) = -2ari(°) = (0,0,fTT,a4)ii;, (77) 

where the spin-resolved Drude conductivity of majority 
and minority spin channels is given as usually by cr|,| — 

However, n^^^ (x) has an additional potential drop that 
is extracted from the asymptotic behavior and that stems 
from the second term in Gii{x), 

voe^SV = nf\x +00) - n^^\x ~oo) (78) 
= 2n'-^\x 00) . 

Here, we used that the charge component ric = ni; + nj^ 
and the second equality is due to symmetry of the con- 
tact. Keeping the current constant, the presence of the 
wall implies a change in the externally applied potential, 
which directly translates into a relative change in resis- 
tance. Hence, we define the DWR 

r- _ Pdw - Po SV 2 ny{x 00) 

Po y ni''\+L/2)-ni°\-L/2) 

(79) 
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Of course, due to the perturbative nature of our treat- 
ment, the correction n^^'> has to be always smaller than 
n'^^K Since it turns out that 

rj(2)(a._,±oo) =(5pDW ^(")(±f) (80) 

this condition is equivalent to (5/3dw 1 which is always 
the case. 



D. Results 



The correction to domain-wall resistance up to order 
= {dm{x))^ takes the form 



6p 



DW 



PDW — PO 
PO 



IE, 



DW 



3 Ef 



/(A,7,^,k) (81) 



where analogously to Ref. [20|, we use the domain wall 
energy defined by 

The constant C of order unity depends on the specific 
form of the wall and we find the geometric scaling to be 
1/wL. For the domain wall profile of eq. we obtain 
C = 2. The scaling with 1/wL can be easily understood 
by realising that corrections to the resistance arising from 
a gradient q yield the behavior Spuw on oc -\. Due 



to physical reasons there cannot be a corrections linear 
in q since the result should not depend on the sign of q, 
i.e. the sense of rotation of m. Since the total length of 
the contact is L and the domain wall constitutes only a 
fraction w/L, the correction for the whole contact should 
indeed be Sg oc l/wL. 

A thorough investigation of the whole hierarchy of 
equations reveals that the result obtained for (5/Odw 
is valid for wall lengths much larger than the spin- 
precession length, w 3> Iprec- The mathematical reason 
for this condition is, that, even though in the vicinity 
of the domain wall, the detailed profile of the quasi- 
particle excitations n^^^x) depends on the whole hier- 
archy of equations (unless w ;» Isd, with the transverse 
spin-diffusion length l"^^ = D±/'y:^), the asymptotic be- 
havior of the correction n^^^(a;), is not affected by higher 
order contributions of the multipole expansion. And in 
situations, in which finite size effects from the contact 
geometry are negligible, Spuw is solely determined by 
these asymptotics. In short, to obtain the asymptotics, 
and thus the desired result, valid up to order g^, we need 
to take into account contributions up to g'--^^ . 

Up to order x^i the final expression for / can be split 
into three parts, 

/(A...f..) = ^'^;^'t^ff^;3',+ <^'^-^' (83) 



A2(l-72)(l + p(A-7)-7A) 
that take the explicit form 



,r , _ 5(2(A-2P) + A^(A-P))-|-A7(l3A(7A-7^-l)-10A+87)-PA7(l9(7A-l-7^A^)-|-9(A^-3)-l-47^(l-|-2A^)) 
f^i^A,^) = -(l-7^)A(l-A^)+2P(l-7A)(l-27A+A^) 

/.(A, 7) = . (84) 



Let us reconsider the assumptions made during the 
derivation of result ([5T|) , x ^ 1 and w ^ Iprcc along with 
(PI). The former can be rewritten as Ig S> ^prcc where 
Is denotes the scattering mean free path. This in fact 
implies that Zprcc is, besides the Fermi-length, the small- 
est length scale in the system. Note that no assumption 
was made on the relation between Ig and the length of 
the domain wall. However, we made the assumption of a 
diffusive contact which implies that L ^ Ig. 

For completeness, let us also specify the leading order 
correction to the longitudinal current component, 

jI^ - (85) 
jf)(x) = eED^^ol^9{A,^,C,n) 
which, in contrast to dpuw, is only valid for w ^ Igd- 



Obviously, due to conservation of charge current, this 

(2) 

correction constitutes a pure spin-current, js {x). The 
x-dependence is inherited trivially from Eq{x) ~ '^^m ■ 

1. Quasiclassical Regime A -C 1 

Let us first investigate the limit A Ep, commonly 
refered to as the quasiclassical regime. However, due to 
the restrictions imposed upon Ip^cc and discussed above, 
A cannot become arbitrarily small. In this regime, cor- 
rections to the electron density of states play no role and 
the density of states can be considered constant. Also, 
it turns out that gradient corrections to the self-energies 
play no role so that only the contribution from impurity 
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scattering remains. We find in this limiting case 

v2 //I \ ]^ 



/(a,7,C,«) = 



7 



1 



■ 7 



(86) 



with the scattering asymmetry parameter 7 as well 
as spin-dephasing strength ^. This result shows the 
strong enhancement by scattering asymmetry as already 
noted in previous works^S,. It also displays the same 
A-dependency already found in various other works. 
Clearly, this result is due to different conductivities in 
the two spin channels which are mixed in regions of non- 
vanishing magnetization gradient q. A spin-up electron 
incident on a domain wall attains a spin-down component 
since electron spin direction does not instantaneously fol- 
low the local magnetization direction m (known as spin 
mistracking) . Since the electronic spectrum plays no role 
in this limit, any asymmetry in the conductivity of the 
two channels is due to the scattering asymmetry 7, thus, 
there is no DWR as 7 0. 

Comparing this result with the work a^^i^°'^^ , we find 
in this limit, besides ^, a different numerical prefactor. 
Adopting to our notation, they have / = | ^Z^2 j?- To 
fathom this discrepancy, we stress that the specific form 
of the longitudinal Greens function Gii(a;) is a result of 
the presence of spin-flip processes. In absence on these 
processes (which is the case in aforementioned works), 
the spin-diffusion length diverges, so that properly per- 
forming this limit yields the longitudinal Greens function 
Gii{x) = He'^ + ^ (K-H)^-^, which produces 
qualitatively different results. This additional contribu- 
tion persists arbitrarily far away from the domain wall 
and leads to a different result for (5/Odw- Calculating 
then the DWR in the quasiclassical limit, we find that it 
is still smaller by a factor of 2 as compared to the result 
iii2£. Nevertheless, this shows that spin-flip processes are 
crucial and cannot be ignored, since the absence of the 
latter leads to spin accumulation that does not decay 
even infinitely far away from the domain wall and thus 
yields an additional contribution to the DWR. 

Concerning the work of^*, the main critic has been 
mentioned in the introduction. Even restricting ourselves 
to the quasiclassical regime, the use of a system consist- 
ing of an infinite spin-spiral the inclusion of only up to 
p-wave component and the lack of terms due to gauge 
transformation can be invoked to explain the discrepancy 
to our fully microscopic results. 

There is still the question about how small A can be- 
come, since, clearly, the limiting factor appearing in the 
expressions is 1/A^. However since, within our specific 
model. 



7 cx P(A) = - -f 0[A]3 , 

there is no problem concerning Spuw, because then A^ 
cancels out. For the current, only A cancels, but there 
is another constriction leading to E'g <^ A which again 
stems from the requirement that Iprcc is smaller than the 










.1 


- ? = 


.4 




A = , 



FIG. 2: Domain- wall resistance 5pDW = Isf' /(^>7i^''^) 
as a function of exchange splitting A for various values of ^ = 
the ratio of spin-dephasing rate and momentum relax- 
ation rate and denotes the transverse spin-dephasing strength. 
E-Dw = ^jj is the domain wall energy, k, — 1 implies spin- 
isotropic magnetic scattering of strength x\\^^ = (dotted 
lines), while k — 2 means Xi™' ~ (normal lines) and k = 

is the case of Xx"' = (dashed lines). 7 = P(l — ^k) is the 
scattering asymmetry parameter and is not an independant 
parameter in our model. 



scale of the magnetization gradient. Hence again, we 
will not get into trouble. This remains a problem in the 
work of Brataas and coworkers^-, since there, the impu- 
rity scattering times are introduced phenomenologically 
and thus, in their scenario no restriction is placed upon 
7- 



2. Arbitrary A < 1 



Let us now have a closer look at the behavior in the 
whole range of valid values for A. In this regime, the gra- 
dient corrections in the collision integral become impor- 
tant and so is the influence of the magnetization gradient 
on the electronic structure, viz., the density of states. 

In the half metallic limit the spin-flip length becomes 
arbitrarily small, since A ^ as A ^ 1. Writing A = 
Ep — and owing to the condition that no length should 
exceed the precession length /pi-oc *C A, we obtain the 



requirement that 



prcc 

» 7,^7,^ 



However, this is not a 



k(1-7) 'n 'J 

big restriction considering our assumption that scattering 
is weak, so that 7^"" , 7;|" Ep. 

Note that within our model of impurity scattering we 
have only two independent parameters. A convenient 
choice is to vary ^ and k independently, so that 7 = 
P(l — ^k) is fixed, k can be also regarded as magnetic 
scattering asymmetry parameter since we have k — for 
= 0, K = 1 for Xii™' = Xj™^ a-nd k = 2 for Xy™'' = 0. 
For various values of ^ and k, the DWR is shown in figure 
[5] and the corresponding change of the spin-current in the 
contact is depicted in Fig. [31 
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s(a,7,C'«) 
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.4 




FIG. 3: Additional magnetization-gradient-induced flow 
of longitudinal spin-current inside the domain wall where 
Eq{x) = \jf^ does not vanish, so that ji^\x) — 



eED±vo 



1 Eg(x) 

3 Bp 



(/(A,7,^,k). Note that the divergence for 



A poses no problem, as discussed above. The magnetic 
scattering parameters are identical to FIG. [2] 



We can see that the DWR does not vary strongly 
with A when this parameter is small and decreases 
monotonously as A increases even to the point where the 
DWR can become negative as one approaches the half 
metallic regime. The latter happens only if Xj^"' 7^ 0, i.e. 
in the case of non-vanishing spin-flip scattering between 
the two bands. Also, the magnetization-gradient-induced 
longitudinal spin-current (c.f. Fig. [3]) displays a qualita- 



tive difference between the cases where 



X± 



(m) 



where it does not. This clearly is a band structure effect 
and we find that this also requires to include corrections 
to the collision integral. It turns out that the presence of 
a magnetization gradient modifies the density of states 
so that there are corrections of order to the scattering 
rates that decrease the momentum relaxation rate and 
a=2If thus, a reduction in resistivity. In fact, the density of 
states for the minority spin channel strongly decreases 
as one approaches A — > 1. Furthermore, we find that 
the neglect of said corrections leads to a monotonic in- 
crease of the DWR with increasing A as opposed to the 
present result. As a side note, we remark that it is im- 
portant to include corrections to the collision integral up 
to the same order as in the approximation of the trans- 
port part of the kinetic equation for the calculation to 
be consistent. Finally, predictions for a given material 
would require realistic band structure calculations, but 
the present calculations demonstrates the possibility to 
have a negative DWR. 

Let us finally compare our approach to the one of 
Brataas et aS2, by temporarily excluding spin-flip pro- 
cesses. In that case, we obtain for the domain wall re- 
sistance a strictly monotonic increase with A compatible 
to the findings of Brataas et al, in particular, there is no 
negative DWR. Another puzzling fact appears, when we 
let the difference in momentum scattering rates vanish, 
by setting 7 = in our result (let us assume for a moment 
that 7 is an independent phenomenological parameter), 
corresponding to r| = in the result of Brataas et aU^. 
vanishes an We find 



/(A,7 = 0,e = 0,K = 0) = 



A^ 



P 4 + A^ 



2A(1 -hPA) 
I 



= A^ 

16 



o[K] 



(87) 



which yields an overall leading order term proportional 
to A^. In contrast the result of Brataas et al yields 
a constant contribution in that limit: 5p^ — | . 
Interestingly, exclusion of spin-flip processes in our ap- 
proach yields a result with the same asymptotics, i.e. 
/ (a, 7 = 0, ^ 0, K 0) = I 0[A]2. Unsurprisingly, 

the results still do not coincide exactly, since we used the 
fully microscopic collision integral in our approach and 
corrections to the electronic structure are already present 
in this order. 



IV. CONCLUSION AND OUTLOOK 

We have calculated the domain-wall resistance of a 
Bloch wall situated in a wire of quasi one-dimensional 
geometry. The current flow is perpendicular to the wall. 



We assumed that any spin-accumulation has decayed at 
the contact ends, thus neglecting any finite size effects. 
Going towards the half-metallic regime, our calculations 
show the existence of a negative DWR in the presence of 
non-magnetic scattering giving rise to spin-flip scatter- 
ing. This possibility to obtain a negative DWR is a band 
structure effect. Disagreement is found when compar- 
ing our results to various previous works on DWR. We 
believe these discrepancies seemingly arise, on the one 
hand, from the neglection of spin-flip processes and, on 
the other hand ,from the different approaches to include 
impurity scattering. While we use a fully microscopic 
approach for scattering, reflected in the full form of the 
collision integral (|39p with gradient corrections and mod- 
ification of the electronic structure properly taken into 
account, other works introduced momentum scattering 
rates phenomenologically. 

To summarize, we have derived fully miscroscopic 
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equation for the spin transport in noncoUinear magneti- 
zation textures. Our approach takes impurity scattering 
and spin-flip scattering into account on the Hamiltonian 
level. This paves the way to treat more complex mag- 
netic textures and derive microscopic expression for the 
domain- wall-induced resistance. 

All previous works dealing with DWR in the limit of 
wide walls obtain results that depend in a similar way 
on the microscopic parameters, that is, the DWR is Co? 
where C is a dimensionless prefactor and a = is 
the spin mistracking angle. Hence, every theory predicts 
this sort of dependency, but with differnt proportionality 
factors, both in value and sign. This factor C however 
contains information about the scattering and DOS of the 
two spin channels and can depend in a complex manner 
on these properties. The most simple one is due to a 
model by Levy and Zhangi^ where C simply depends on 
the ratio of resistivities for spin up and down channels. 

We did not take into account the possibility of mag- 
netic moment softening, i.e. the reduction of magnetic 
moment within the domain wall. This effect is most 
prominent in very sharp domain walls where canting of 
adjacent spin is large so that the noncoUinear spin states 
hybridize which in turn leads to a reduction in the ab- 
solute value of the magnetic moment. As shown in^, a 
reduction of the magnetic moment can lead to a negative 
DWR. 

Finally, effects due to geometric confinement have not 



been considered, for example, surface scattering might 
become important. Also, the magnetization profile can 
be more complicated and might lead to eddy currents in 
the vicinity of the domain wall which might be relevant 
for interpretation of experimental results on the DWR in 
thin nanowircs. 
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APPENDIX A: APPENDIX - DETAILS ON 
DERIVING THE HIERARCHY OF EQUATIONS 

We start out with the kinetic equation ([55)1 and 
by taking (— e) (u" p8|) ) and by introducing 

a"" = 9) = / 9 , 

Afterwards, employing an expansion of o = 
ei(axOfc,-afe, a J ~ 1 i(- • • ) - i(- • • )2 and making 
use of relation (|49l) , we obtain 



fix) ° g{x, k) 



f,9 



(n-l) 



8A/2 2 
n(n — 1)1 



?i("-2) 



SAP 



(Al) 
(A2) 



where we are able to truncate the series since our treatment includes only terms up to order and subsequent terms 
would contribute only to higher orders. 
This yields the following equation 



dg 



(«+i) 



m(T,g 
1 



(") 



An 1 
2M2 



dm 



■'} 



An{n — 1) i 



8M 



42 



d^ma, g 



(n-2) 



i({a("-^),92#(")}-{a^iaW,3("-^)} 



n{n — 1)1 

where the action of x is defined as 

±X = {x^^ + Trx^™)) [X -iKX\\a-^ x^a 



X — S^K x\\ar — x^a 



r 



and we write x\\ = Tr{mo-A"}m and a;_L = Tr{o-A"} — a;|| , simply becomes 

so that the total impurity self-energy, consisting of spin- 

isotropic and magnetic parts, equations (jl9[) and (|20p . -y _ 



Si 



(27r)3 ^ 



(A3) 



(A4) 
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We also need to know the spectral density a*^"^ up to order q^. We find for even and odd moments, respectively 
I 



^(2«+l) 



27rQ;||"-*m X 9m(T 



(A5) 



where {d^m)± denotes that we only take the component i/^^i specified without any argument implies that we take 



perpendicular to m. The coefficients for a d — 3 dimen- 
sional electron gas are given by 



An) 



A 



J_ („) _ 2n(2n - 1) („_i) 



8Af2 



2n-l („_i) 2n-l fn-i)6fi 
a I H i^k ^ 



1 (n) 

'2A"-^ 



2n - 1 
8M ' 



1 - 



9 

8M 



its value at the Fermi-level, i.e. more specifically v 



using the definitions 



1 



T^J- ^ ^ 1 + 2n 



2 / A 



The chemical potential Sfi is obtained from the condi- 
tion that 



+ C30 



^/DHTr{a(°)H-a(°)H|,=o} 



+ 00 



+ 00 



M \4 

since in the zero-temperature approximation 



— OO 

1 

— ( 

2 



(A6) 



4-00 /"-Ef 



M 



(n+l) 



2n + 1 



Therefore, we immediately arrive at 
5fji 



1 



2 (0) _ 1 2 Af4'^^ 



(0) -L 2Af 4 



- 1 



A 



, (A7) 



which yields equation pop . once we plug in all defini- 



tions and by noting that 



(0) 



(0) 



Pvq and 



,(1) 



■^{2PEy + A). Substituting back into the 



expressions for /3, we can now write 



2n- 1 
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(n-l) 
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Later, we will also need 

iaW = ±a^°^+±a^°^ ^ + P{1 - C«)m^) + ^ i + /3i°)(l - ^n)m&) q' + - 0{d^m)^& , 



where we write ciq'^ to indicate that we take the zeroth or- 
der in q only and accordingly, Uq^^ contains all corrections 
to the density of states due to a magnetization gradient. 

Next, wc will change the 2x2 spin-matrix representa- 
tion to the 4x4 matrix representation in the basis intro- 
duced above, {+, — , T, i)- This change of basis implies a 
local rotation of the basis in spin space to align the mag- 
netization direction m along the new z-axis, which cor- 
responds to the gauge transformation introduced above. 
The various vectors we will encounter in the following 
derivation transform in the following manner: 



and we note that —md^m = {dm)'^ = q^- In the fol- 
lowing, we specify substitution rules that perform this 
transformation: 



ai,n 



\ {ai,n} ^ Kj,n 



(A8) 



ma 

dmcr 

m X dm 

d^mcr = {d^m)±a + {md^m)m& 



qay 

-qo-x 

{dq)ay 



where 




K„ = 



/ 0-11 
0-11 

110 

V -1 -1 

1 1 \ 



110 0/ 




-z \ 
z 

\ / 
\ 

10 
-1 / 



constitute matrices in (+,— ,T,i) representation. The As a consequence of the gauge transformation, the 
action of x turns into derivative transforms into 



7, 



Xn 



with the matrix 



d^d=ld + qM^ 



-10 

0-10 
-f f 

f -f 



Now with these rules at hand, the change of represen- 
tation is straightforward and we obtain 



■2) 



(A9) 



where zeroth order relaxation and precession terms are 
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and magnetization gradient correction to relaxation rates yield 

Corrections that depend on lower moments in the hierarchy read explicitly 

An 1 



I {ia(°\.9(")} - ^ [g^/3^"^l + - ^^)K. + (9<z)/3f (1 - ^Ky 



2A/ 2 I ' ^ i Ar ^'^ 2M 2 L - ' ^ 

(5<5M)l + g:^7^7M. 



o if,, 



n 



and 
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n(n — 1) 
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8M 



2 [A - idq)My) - 7/7 - idq)^y 
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8M2 2 
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(All) 



(A12) 



(A13) 



Furthermore, we have source terms appearing in the for even indices and ^q^"^^^ = for odd indices. The 

corresponding gradient corrections are, for even indices 



equation, whereof the zeroth order term simply is 



^^^0) ^ .(0) 

(A14) 
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8A'P 
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and for odd indices 



1 //j(2n+l) 
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^ / 2M 2 
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^ ^ ^ 2M " = = 



^0) ^g(2n+l)^0) ^ (A16) 
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For our purpose, we only need the first 5 equations, and U = g^^\ Here, we already dropped terms that 

since, as stated previously, in our regime of investigation would only contribute to higher orders than q^. Note 

we need to know only up to g^"*'. Explicitly, these equa- that F = 11 — Sq*'-'. 

tions read Our aim is to obtain a differential equation of the form 





C"(o) _ 


n("))n 




(A17) 


dS + Uj = 
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df + nS-E^'>n = 




-nf 5- 




(A19) 


dU + Uf = 






-0{q') 


(A20) 




0{q) , 






(A21) 



(F - Dd^) n = W{q)n 



(A22) 



where 



W{q)=Y,oqd + Y,,dq+ 

Y^o,,adq^d + Y^,a„oqdqd + Y^o,aqdqdq. (A23) 

To achieve this, we unite the set of equations (jA17P - (|A2ip 
where we defined n — (f^\ j = g^^\ S — g^^\ T = g'^^^ iteratively by eliminating every moment except n. We 
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do this order by order in q,d and in the following, we 
give only terms that are relevant to our result. The ze- 
roth order is simply Sq = H^^Sq^^ti and Uo = H^^Sp'^^n 
while the odd moments j , T vanish as q, 9 — > 0. Sq, when 
plugged into equation (|A18p will yield a term that resem- 
bles a diffusion term, but will not yet be of the desired 
form as given by (j55p . In order to tranform the diffu- 
sion term into the form which permits convenient solu- 
tion of the differential equation, we will multiply equation 
(|A22[) to the left with HA and then subtract it from the 
equation (|A19p which provides S. We thus obtain a new 
•So — (^n^-^sj,^^ + n with the choice of A determined 

by the condition that = WD n which in turn translates 
into 

A = (n£ - n-^S^^') (A24) 

and, since F is singular in longitudinal subspace as it pos- 
sesses one vanishing eigenvalue, F^^ has to be considered 
as the pseudo inverse. 

In order to obtain the form reminiscent of the spin- 

I 

The first order is given by (j = {J_qq + J_ad) n, etc.) 



charge diffusion equation (|50p . we rewrite equation (jA18|) 

U{j+Ddnj = UDdn - 95 + Sj^^n - n9-'>i+&^n , 

(A25) 

where now the right-hand-side of this equation is of or- 
der q, in particular the difference HP n — S no longer 
contributes to the zeroth order solution, as opposed to S 
itself. 

Before proceeding with higher orders, we define for 
convenience 



n(") 

=9 
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kl 
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= dtl'+dtidq-dq) 


2(2") 


- m'''<i'+^Vdq+^f 


2(2"+l) 





where we remind that d acts on everthing to its right so 
that actually are differential opperators. 



the second order terms are 



and finally, the 3 remaining coefficients in equation (|A23p . 
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qrlqrl 
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To solve equation (jA22p . we use the method elaborated upon previously in this article, so that by use of equations 
74]) and we end up with 



<5pDW 



. [2(rajit£(r,a)tl + (^a,)lt£(^ajtl + i2Ya,,a +Y 

qdqd + ^3,39)11] '^0 



(A26) 
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where the outer bracket denotes that we take the charge co mponent. 
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